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I.  INTRODUCTION 


Flaxe-stabilized  projectiles  are  currently  being  evaluated  for  a  variety  of  applications, 
including  a  25  mm  training  roundd  A  systematic  study  of  a  number  of  flare  geometries 
has  been  undertaken  at  the  U.S.  Army  Ballistic  Research  laboratory.  The  ultimate  objec¬ 
tives  of  the  study  axe:  1.)  to  determine  how  to  increase  the  stability  of  a  flare-stabilized 
round  with  a  minimrmi  increase  in  drag;  and  2.)  to  provide  a  range  of  aerodynamic  data 
with  which  to  determine  the  capabilities  and  limitations  of  available  numerical  predictive 
techniques. 


1.  DESCRIPTION  OF  THE  CONFIGURATIONS 


Ten  projectile  configurations  are  considered  in  this  investigation,  all  with  basically 
the  same  slightly  blunted  cone- cylinder  forebody  of  almost  12  calibers.  The  projectile 
afterbody  consists  of  a  stabilizing  flare,  the  angle  and  length  of  which  were  varied  to 
provide  a  range  of  flare  shapes.  Figure  1  shows  the  overall  geometry  of  the  projectile  with 
one  of  the  flares  (CS-V4-1)  attached. 

Figure  2  shows  all  of  the  flare  geometries.  The  CS-Vl  configuration  is  a  simple  15.2 
degree  half-angle,  2.67  caliber  long  flare.  The  CS-V2  design  has  a  4  degree,  4.5  caliber 
fiax'e.  The  CS-V4  family  of  flares  are  all  basically  derived  from  a  6  degree,  3.5  caliber 
baseline  flare  configuration.  Configuration  CS-V4-1  is  the  baseline  flare.  Configurations 
CS-V4-2  through  CS-V4-5  consist  of  the  baseline  flare  with  one  caliber  extensions  of  6,  0, 
12,  and  -6  degrees,  respectively. 

Additionally,  there  a-e  three  other  configurations  in  the  CS-V4  group  that  are  not 
simple  conic  or  biconic  shapes.  CS-V4-7  has  a  -6  degree  boattail,  similar  to  CS-V4-5, 
but  it  also  has  a  set  of  four  fins  with  a  12  degree  sweep  angle  located  on  the  boattail. 
Configuration  CS-V4-8  consists  of  the  baseline  configuration  with  four  strakes  extending 
along  the  length  of  the  flare. 


Finally,  configuration  CS-V4-6  has  a  non-conical  square  base.  The  geometry  of  this 
shape  was  chosen  so  that  the  base  area  would  be  identical  to  that  of  the  CS-V4  baseline 
and  so  that  the  last  caliber  of  length  would  have  a  square  cross  section.  This  dictates  a 
flare  angle  of  9.53  degrees. 


Note  that  the  basic  shape  corresponds  to  that  of  anti-armor  projectiles  of  relatively 
large  length  to  diameter  ratio  (about  16.3)  where  the  center  of  gravity  is  near  the  mid¬ 
length.  In  order  to  maintain  the  same  weight  distribution,  the  projectile  bodies  were  made 
of  steel  and  the  afterbodies  were  made  of  aluminum  with  an  open  cavity  in  the  base. 

The  designations  of  the  ten  flares  and  flare-fin  stabilizers  and  their  major  geometrical 
relationships  are  given  in  Table  1. 
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2.  EXPERIMENT 


A  series  of  test  firings  have  been  carried  out  in  the  BRL  Aerodynamics  Range. ^  Up 
to  four  rounds  of  the  CS-Vl  and  -V2  projectiles  were  fired  at  a  number  of  fixed  conditions 
from  transonic  to  4.6  in  Mach  number.  Four  rounds  of  each  of  the  other  flare  configurations 
were  fired  at  a  nominal  Mach  number  of  4.0.  The  sabots  used  to  launch  the  projectiles 
incorporated  a  slip  barid  obturator,  which  is  designed  to  keep  the  roll  rate  at  approximately 
one  turn  in  four  meters  of  travei.  Orthogonal  spark  shadowgraph  pictures  (see  Figures  3 
and  4  for  typical  pictures)  at  stations  along  the  range  are  interpreted  to  provide  position 
and  orientation  throughout  the  flight.  These  data  are  fitted  to  the  linearized  equations 
of  motion.3  fitting  provides  the  linearized  aerodynamic  coefficient  of  drag  and  the 
slopes  of  the  coefficient  of  moment  and  lift  versus  angle  of  attack.  From  the  lift,  the 
center-of-pressurc  is  also  determined. 

TABLE  1:  GEOMETRICAL  PARAMETERS  OF  THE  XM910  CS-FAMILY  OF  PROJECTILES 


CONFIGURATION 

FLARE 

FINENESS 

RATIO 

FLARE 

ANGLES 

BASE 

DIAMETER 

RATIO 

MAX. 

SPAN 

RATIO 

EXTENSION 

ONE 

CALIBER 

CS-Vl 

2.67 

2.450 

None 

CS-V2 

4.50 

4.00 

1.628 

None 

CS-V4-1 

3.50 

1.736 

None 

CS-V4-2 

4.50 

1.948 

6.0°  Flare 

CS-V4-3 

4.50 

1.736 

Cyl.  Skirt 

CS-V4-4 

4.50 

2.160 

12°  Biconic 

CS-V4-5 

4.50 

1.526 

-6°  Boattail 

CS-V4-6 

4.50 

9.53 

1.540 

Square  Base 

CS-V4-7 

4.50 

1.226 

2.160 

Boattail  w/  Fins 

CS-V4-8 

3.50 

6.00 

1.736 

2.043 

None 

II.  DESIGN  CODE  ANALYSIS 


Two  kinds  of  predictive  techniques  have  been  applied  to  the  family  of  configurations: 
feist  turnaround  design  codes  which  are  ideally  suited  to  parametric  studies;  and  the  Parab¬ 
olized  Navier-Stokes  (PNS)  technique  for  solving  the  supersonic,  thin  sheai'  layer  equations 
of  motion. 


1.  SELECTION  OF  A  DESIGN  CODE 

The  Missile  DATCOM**  and  the  NSWC  Aeioprediction^  codes  are  cap.ible  of  analyzing 
the  aerodynamic  characteristics  of  the  proposed  XM910  configurations.  The  two  codes  use 
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the  same  two  methods  for  the  evaluation  of  the  inviscid  contribution  to  the  aerodynamic 
coefficients  at  supersonic  speeds.  They  are  Van  Dyke’s  hybrid  method*'  and  Syvertson 
and  Dennis  Second  Order  Shock  Expansion  (SOS£)  method.'  Indeed,  both  codes  yield 
practically  the  same  results  for  the  normal-force  and  pitching-moment  coefficients  at  small 
angles-of-attack  and  the  same  zero  lift  wave  drag  coefficients. 

In  the  transonic  region,  the  codes  use  data  bases  for  the  evsduation  of  the  normal-force 
and  the  pitching  moment  coefficients.  The  accuracy  of  the  two  data  bases  weis  evaluated 
by  comparing  predictions  by  the  two  codes  with  the  experimental  data  of  Barth.®  The 
configuration  selected  for  the  comparison  is  a  cone-cylinder  having  com,ponent  fineness 
ratios  of  2.5  and  6.0,  respectively.  This  coiiliguration  is  similar  to  that  of  the  forebody  of 
the  subject  round. 

The  results  of  the  comparison  are  shown  in  Figures  5  and  6.  The  slope  of  the  normal- 
force  curve  predicted  by  the  DATCOM  code  is  in  excellent  agreement  with  the  experimental 
data  base.  The  NSWC  code  predicts  a  slope  which  is  11  to  27  per  cent  lower  than  the 
experimental  values.  The  DATCOM  code  predicts  the  location  of  the  center-of-pressure 
within  0.25  diameters  and  better  than  the  NSWC  code.  These  findings  concerning  the 
data  base  of  the  NSWC  code  sire  consistent  with  previous  conclusions  of  Sigal.®  It  wais 
decided  to  use  the  DATCOM  code  for  the  evaluation  of  the  aerodynamic  coefficients  in 
the  pitch  plane. 

The  two  codes  use  different  empirical  methodologies  for  the  estimation  of  the  viscous 
contributions  to  the  axial  force  coefficient.  Since  the  outputs  are  different,  these  contri¬ 
butions  were  compared  with  the  empirical,  data  reported  by  Stoney.^'^It  was  found  that 
the  DATCOM  code  ove/-predicts  the  base  drag  coefficient  of  bodies  having  flares.  The 
NSWC-Aeroprediction  code  provides  reliable  prediction  of  the  base  drag  coefficient  for 
cylindrical  and  flared  bodies.  It  uses  an  empirical  data  base  for  the  estimation  of  the  skin- 
friction  coefficient  in  a  wide  range  of  Reynolds  numbers,  including  the  region  of  transition. 
This  capability  is  of  importance  in  the  present  case  since  the  small  projectile  flies,  most  of 
the  time,  in  the  transition  zone.  The  NSWC  code  was  selected  for  the  evaluation  of  the 
axial-force  coefficient. 


2.  CHOICE  OF  DESIGN  OPTIONS 

In  the  beginning,  the  two  options  for  calculations  at  supersonic  speeds  were  tried,  in 
em  attempt  to  find  the  best  cross-over  Mach  number  between  the  Van  Dyke  and  the  30SE 
methods.  However,  when  the  results  were  analyzed,  it  was  found  that  the  Van  Dyke  method 
failed  in  the  case  of  bodies  with  flares.  The  computed  pressure  distribution  along  the  flare 
showed  a  dip,  followed  by  a  rise.  This  distribution  is  physically  unreasonable.  Hence,  only 
the  SOSE  method  was  used.  Since  this  method  is  valid  mainly  for  high  Mach  numbers, 
it  is  desirable  to  evaluate  it’s  accuracy  at  the  low  end  of  the  supersonic  region.  This  was 
done  by  comparing  analytical  results  for  an  extended  nose  lorebody  with  empirical  data 
bases.  The  British  Engineering  Science  Data  Units  (ESDU),“  Aircraft  DATCOM, ^^ajid 
AMICOM,^®data  bases  were  used  and  the  results  are  shown  in  Figures  7  and  8.  The  first 
figure  shows  that  the  SOSE  method  is  in  very  good  agreement  in  normal-force  cuive  slojie 
at  Mach  numbers  larger  them  3.0  with  the  ESDU  and  Aircraft  DATCOM  data.  In  the  lower 
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range,  the  method  over-predicts  the  empirical  data  b}'  about  20  per  cent.  The  center-of- 
pressure  location  is  well  predicted  by  the  SOSE  method.  Note  that  since  the  AMICOM 
data  is  only  a'vai’'blc  for  the  nose  fineness  ratio  of  4.0,  it  was  adjusted  by  multiplying  the 
center-of- pressure  distjuice,  measured  from  the  tip,  by  the  ratio  of  nose  fineness  ratios. 

An  attempt  to  use  the  truncated  nose  option  in  the  code  yielded  a  very  large  reduction 
in  normal-force  curve  slope,  accompanied  by  a  forward  shift  in  the  ceiiter-of-pressure. 
These  changes  axe  a  result  of  the  code  accounting  for  losses  through  a  normal  shock  \'.  ave 
and  not  considering  the  mixing  in  the  entropy  layer.  It  was  decided  to  extend  the  nose,  in 
the  computational  model,  to  form  a  complete  cone  and  to  apply  corrections  for  the  effect 
of  bluntness,  based  on  empirical  data,  that  will  be  described  in  section  5. 


3.  PRESENTATION  OF  THE  NON-CIRCULAR  STABILIZERS 


Three  configurations  feature  non-circular  stabilizers.  Configuration  CS-V4-0  has  a 
square  base  and  CS-V4-7  and  -8  have  lifting  .surfaces  mounted  on  the  flares.  Neither  feature 
can  be  handled  by  the  design  codes  presently  availaole  at  the  Launch  and  Flight  Division. 
Hence,  equivalent  axisymmetric  configurations  were  used  as  computational  models  in  these 
crises.  Their  cross-sectional  area  distribution  provide  the  same  normal-force  curve  slope 
distribution  as  that  of  the  actual  configurations,  evaluated  by  slender  body  theory. 

Using  center  body  cross  sectional  area  as  reference  area  for  the  equi\^Jent  body,  the 
above  condition  becomes: 

(^).  =  ^  (1) 


5/z 


Sr 


Since  the  equivalent  body  is  circular 


^  Sr 


(2) 


For  a  slender  body  having  a  square  cross-section,  of  dimension  d  on  a  side,  Nielsen, i-igives: 

5Cw<,  2Ad^  _ d 


Sr 


Hence,  for  a  square  afterbody 


-^-  =  3.056(~)^ 


For  a  cruciform  slender  wing-body  combination  the  same  reference  gives 


so  that 


D  o'-  ^  b  ^  ^  b  ^  ^ 


(4) 


(5) 


(6) 


The  main  geoinetricad  parameters  for  the  equivalent  afterbodies  are: 
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CONFIGURATION 

DJD 

x/d-=15.65 

X/D=1G.C5 

CS-V4-6 

■HH 

CS-V4-7 

CS-V4-8 

A  graphical  comparison  between  the  equivalent  and  the  actual  cifterbodies  is  shown  in 
Figure  9. 

4.  THE  DR  AG  OF  THE  BUTTRESS  GROOVES 

This  contribution  to  the  axial-force  coefficient  is  not  included  in  the  design  codes. 
It  is  evaluated  using  the  method  of  equivalent  sand  roughness  and  a  new  correlation  for 
two  dimensional,  ridge-type,  roughness  elements,  compiled  by  Sigal.^^The  ratio  of  the 
equivalent  sand  roughness  to  the  height  of  the  buttress  grooves  is  estimated  to  be  2.0.  The 
amplification  in  the  skin-friction  due  to  the  roughness  is  evaluated  using  a  friction  law  also 
compiled  in  Reference  15.  The  additional  drag  coefficient,  due  to  the  grooves,  is  related  to 
the  amplified  skin  friction  by: 

=  (7) 

i>R  (-'fo 

The  dependence  of  this  contribution  upon  Mach  number  depends  on  the  wall  temperature 
conditions.  A  cold  wall  case  is  assumed  here,  where  the  wall  temperature  equals  the  free- 
stre£un  temperature  because  it  approximately  represents  the  conditions  in  the  aerodynamic 
range  tests,  where  the  time  of  flight  is  too  short  for  pronounced  aerodynamic  heating.  In 
the  transonic  region  the  difference  between  a  cold  and  adiabat  .  wall  conditions  is  small. 
As  Mach  number  increases  the  adiabatic  wall  case  yield  a  smaller  contribution,  due  to  the 
increase  in  wall  temperature. 


5,  CORRECTION  FOR  THE  TRUNCATED  NOSE 

As  discussed  above,  the  nose  in  the  computational  model  was  extended  to  form  a 
complete  cone.  The  correction  due  to  the  actual  geometry  is  based  on  empirical  data. 

The  systematic  data  base  of  Owens, ^^‘for  spherically  blunted  cones,  is  used  for  the 
correction  of  the  normal-force  and  the  center-of-pressure  location.  Results  for  a  half  cone 
angle  of  10  degrees  are  used  because  they  are  closest  to  the  present  geometry.  The  effects 
of  bluntness  is  a  slight  decrease  in  normal-force  curve  slope  and  a  very  slight  forward  shift 
in  the  center-of-pressure. 

Since  no  data  were  found  for  truncated  cones,  the  estimate  of  the  effect  of  truncation 
on  the  axial-force  coefficient  is  based  on  the  data  of  Escli^'  for  truncated  ogives. 

The  equivalent  configuration  used  for  the  evaluation  of  the  contribution  of  the  skin- 
friction  and  the  base  to  the  drag  coefficient  have  the  same  wetted  area  and  base  area,  as 
the  actual  ones. 
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In  the  transonic  region,  the  design  codes  use  slender  body  theory  for  the  rvahiation  of 
the  contribution  of  afterbodies  to  stability.  Thus,  only  base  diameter  ratio  (not  the  base 
area  distribution)  determine  the  contribution  of  the  stabilizer  to  the  iioriual-force  curve 
slope.  As  a  result,  the  design  codes  yield  the  same  normal-force  curve  sloiio  and  almost 
the  same  center-of-  pressure  location  for  configurations  CS-V4-1  and  -3,  in  the  transonic 
region. 


6.  COMMENTS 

Configuration  CS'V4-5  has  a  boattail  mounted  on  a  flare.  The  flow  deflection  between 
these  two  components  is  12  degrees.  It  is  expected  that  this  large  deflection  will  cause 
separation  at  transonic  and  low  supersonic  Mach  numbers.  Hence,  it  is  expected  that 
the  boattail  will  be  aerodynamically  ineffective  in  these  regions  and  that  the  aerodynamic 
chajacteristics  of  this  configuration  will  be  close  to  those  of  CS-V4-1. 

The  fin  of  configuration  CS-V4-7  is  mounted  on  a  boattail  identical  to  that  of  CS- 
V4-5.  By  the  same  argument,  it  is  expected  that  the  fin  will  be  immersed  in  separated 
flow.  This  will  reduce  it’s  efficiency  as  a  stabilizer  so  the  aerodynamic  characteristics,  at 
transonic  speeds,  will  be  close  to  those  of  CS-V4-1. 


Ill-  PARABOLIZED  NAVIEn-SiuKiiii?  cJUDJfi 


The  flow  over  these  configurations  also  have  been  computed  using  tl)e  Parabolized 
Navier-Stokes  technique  based  on  the  code  developed  by  Schiff  and  Steger.’®This  technique 
consists  of  generalized  body  fitted  coordinates,  thin  shear  layer,  linearized  ,  approximately 
factored  equations.  A  finite  difference  algorithm  is  solved  implicitly  using  a  block  diagonal 
method  as  formulated  by  Beam  and  Warming. ^^he  approximation,  which  permits  march¬ 
ing  in  the  primary  flow  direction,  limits  the  application  to  attached,  supersonic  free-stream 
flows.  Thus,  the  results  from  the  PN3  code  cannot  provide  information  on  the  separated 
flow  in  the  near  wake  or  the  base  pressure. 

The  code  uses  the  Baldwin-Lomax^i  algebraic  turbulence  model  to  compute  skin 
friction  forces.  The  flow  is  assumed  to  be  turbulent  from  the  nose  although  this  may 
somcwu.at  overestimate  the  actual  viscous  contribution.  A  new  iechnique,’^for  coiuputing 
buttress  tlutad  and  groove  effects  was  applied  to  this  configuration.  However,  the  body 
is  so  small  (S.28  mm)  and  the  flight  Reynolds  number  is  so  high  that  the  technique  is 
not  valid  since  the  estimated  equivalent  roughness  height  is  larger  than  the  undisturbed 
boundary  laj'^er  thickness.  Thus,  the  same  estimate  of  buttress  groove  drag  which  was  used 
for  the  design  codes  was  also  applied  to  the  PNS  results. 

In  principle,  the  PNS  code  can  handle  all  tlie  non-couical  shajK's  considered,  including 
the  fin  and  strake  bodies  hut  these  shapc.s  have  not  been  attempted  here.  PNS  results 
are  reportea  only  for  the  axisyrnmetric  shapes  and  the  square  base  configuration.  I'iie 
tabulated  PNS  results  for  the  wave  and  friction  drag  at  four  Mach  numbers  are  given  in 
Table  2. 
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Some  instability  was  observed  in  the  PNS  caJculations  at  the  cylinder-flare  junction, 
peirticularly  for  the  lovr  Mach  numbers.  This  problem  was  so  severe  for  the  CS-Vl  con¬ 
figuration  with  the  15.2  degree  flare  that  the  results  were  considered  unreliable  and  are 
not  included  here.  The  oscillations  for  the  small  flare  angles  was  reduced  or  eliminated 
by  adding  a  small  fairing  between  the  cylinder  and  the  flare.  In  all  the  results  cited  here, 
a  0.25  caliber  sine  function  fairing  was  u.sed  to  provide  a  continuous  body  radius  and  its 
first  derivative  in  the  trrmsition  region.  Computations  at  Mach  number  4.0  showed  that 
no  fairing  was  required. 


1.  BASE  PRESSURE  AND  BASE  DRAG 

Estimates  of  the  base  drag  for  these  projectiles  were  obtained  from  the  NSWC  code 
which  calculates  the  ba^e  pressure  from  empirical  formulas  which  depend  on  the  equivalent 
body  effective  flare  angle.  In  the  present  case  there  is  some  uncertainty  in  the  application 
to  the  two  element  afterbodies  of  the  CS-V4  family  of  projectiles.  An  attempt  was  made 
to  investigate  this  question  by  using  the  results  from  the  more  elaborate  but  still  semi- 
empirical  work  of  Mueller,  et  al.^°The  accuracy  of  Mueller’s  method  for  flared  afterbodies 
has  not  been  generally  established;  thus,  simple  correlation  formulas  have  been  developed 
to  follow  the  trend  of  his  results  and,  as  a  consequence,  these  predictions  should  only  be 
considered  as  primarily  illustrating  trends. 

The  plotted  results  of  Reference  20  3ugge.st  that  the  base  pressure  can  be  described 
by  the  addition  of  two  terms.  The  first  term  represents  the  base  pressure  of  a  long  (nose 
effects  neglected),  cylindrical  afterbody  which  is  primarily  a  function  of  Mach  number.  A 
simple  power  law  representation  for  this  term  is: 

(  — )cv/.n*r  =  1.9AC-''®  (8) 

Poo 

This  expression  agrees  with  experimental  data  in  the  Mach  number  range  from  3.0  to  4.0 
but  predicts  too  high  a  base  pressure  at  Mach  number  of  2.0. 

The  second  term  accounts  for  the  angle  of  the  flare,  (f>f  and  the  length  of  the  after¬ 
body,  (I//D).  An  cxponsntisLl  expression  is  fitted  to  ^'/Inellei’^s  results.  Tlie  expoiienti3,l 
form  is  chosen  so  that  near  <b]  and  L/  of  zero,  the  approximation  is  linear  and  at  larger 
values  of  these  variables,  negative  pressures  are  avoided  (although  these  formulas  are  only 
approximately  valid  for  small  flare  angles  less  than,  15  degrees,  and  for  Mach  numbers  in 
the  range  from  2.0  to  4.0).  The  additive  effect  of  flares  is  ; 

A(-— )  =  -[1  -  (9) 

Poo 

Note  that  gi/  is  in  degrees  in  this  formula. 

For  the  3.5  and  4.5  caliber  flares,  the  flare  length  term  is  essentially  unity.  Positive 
b/’s  correspond  to  flares  which  reduce  the  base  pressure  below  that  of  a  cylindrical  body. 
Negative  <?/’s  ,  on  the  other  hand  are  associated  with  boattails  which  have  increased  base 
pressures. 

Most  of  the  CS-V4  senes  hav'e  two  element  afterbodies  consisting  of  3.5  caliber,  C 


degree  flaxes,  followed  by  one  caliber  extensions.  The  above  analysis  was  applied  assuming 
superposition  of  the  flare  and  extension  to  the  basic  cylindrical  base  pressure.  However, 
there  is  no  theoretical  or  experimental  evidence  to  justify  this  procedux'e.  The  base  pressure 
is  estimated  as:  n  n 

—  (  )cyl  +  )6“//are  +  A( -^jextenjton  (10) 

Poo  Poo  Poo  Poo 

where  equation  9  is  used  to  evaluate  both  the  6  degree  flare  and  the  extension  contribution. 
In  the  case  of  the  extension,  the  local  body  angle  was  measured  relative  to  the  underlying 
6  degree  flare.  The  fl2ire  length  term  is  essentially  unity  for  the  6  degree  flare  but  does 
have  some  effect  on  the  extension  contribtion. 

Table  2  tabulates  representative  results  based  on  the  PNS  computation  of  the  wave 
and  friction  drag.  The  base  drag  coefficient  is  computed  based  on  the  above  described 
formulcis,  assuming  a  uniform  base  pressure  acting  on  the  physical  base  area. 


IV.  RESULTS 

The  design  codes  (NSWC  and  DATCOM)  and  the  PNS  codes  have  been  applied  to 
the  configurations  of  the  XM910  prototypes  and  the  main  results  will  be  described  in  terms 
of  the  drag,  C^o,  slope  of  the  normal  force  relative  to  the  angle  of  attack,  and  center 
of  pressure,  Xcp.  This  will  be  followed  by  some  results  on  the  detailed  pressure  and  skin 
friction  distributions. 


1.  DRAG  BREAK-DOWN  AND  MACH  NUMBER  DEPENDENCE 


Figure  10  shows  the  variation  of  the  zero  yaw  drag  of  the  CS-Vl  as  a  function  of  Mach 
number.  The  experimental  data  eire  in  good  agreement  with  the  NSWC  codes  prediction 
at  supersonic  speeds  but  the  code  underpredicts  the  drag  at  Mach  numbers  near  and  below 
one.  It  was  not  possible  to  reliably  calculate  this  configuration  with  the  PNS  code  because 
the  15.2  degree  flare  caused  large  oscillations  in  the  calculations.  These  oscillations  are 
believed  to  be  related  to  local  flow  separation.  The  sharp  jump  in  the  predicted  drag 
at  transonic  speeds  (1,15  <  Moo  <  1.20)  is  apparently  caused  by  the  switching  from  the 
low  speed,  data  base,  method  of  computing  drag  to  the  Second  Order  Shock  Expansion 
(SOSE)  method. 


Figure  11  shows  the  Mach  number  behavior  of  the  smallest  flare  angle  configuration 
(CS-V2)  as  well  as  the  components  of  the  drag.  Both  the  modified  design  code  and  the 
PNS  code  predict  the  data  within  the  scatter  in  the  measurements  above  Mach  3.  Belov/ 
that  point  the  scatter  in  the  data  becomes  large  but  with  the  NSWC  again  underpredicting 
the  data.  The  base  drag,  at  speeds  below  about  Mach  2.,  is  equal  or  greater  than  all  the 
otner  contributors.  The  other  factors  remain  surprisingly  constant  over  the  Mach  number 
range  considered  here.  The  PNS  code  results  are  also  indicated  on  the  figure.  The  PNS 
code  can  only  compute  the  pressure  and  friction  components  of  the  drag.  The  base  drag 
has  been  obtained  from  the  correlations  based  on  Mueller’s  technique.  These  calculations 
have  been  added  to  the  PNS  results  as  well  as  the  buttress  groove  drag.  These  composite 
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PNS  results  are  in  better  agreement  with  the  data  than  those  of  the  design  code.  However, 
this  is  in  part  due  to  a  liighep  skin  friction  drag  predicted  by  the  Navier-Stokes  code.  One 
reason  for  the  larger  friction  drag  is  that  turbulent  flow  was  assumed  from  the  steirt  of  the 
numerical  computation  whereas  the  design  code  incorporated  an  empirical  transition  and, 
for  this  small  shell,  this  meikes  a  significant  difference. 


2.  SLOPE  OF  THE  NORMAL  FORCE  AND  CENTER  OF  PRESSURE 

Figure  12  shows  the  Mach  number  dependence  of  the  normal  force.  The  experiments 
indicate,  and  both  computational  methods  predict,  that  there  is  very  little  variation  of 
the  lift  with  supersonic  Mach  number.  For  the  CS-V2  case,  both  the  PNS  and  DATCOM 
codes  predict  the  data  with  equal  accuracy.  The  tendency  for  the  CS-Vl  experimental 
to  increase,  above  Moo  =4.0,  may  not  be  valid  because  of  the  large  scatter  in  those 
particular  measurements. 

The  center  of  pressure  results  are  given  in  Figure  13.  It  is  seen  that  for  the  CS-Vl 
case  the  DATCOM  code  does  well  except  for  a  small  overprediction  at  subsonic  speeds. 
However,  both  codes  tend  to  predict  a  more  forward  Xcp  (by  0.5  to  1.0  caliber)  for  the  4 
degree  flare  ceise  compared  to  the  data.  Thus  the  lift  produced  by  the  flare  is  higher  than 
estimated  by  these  codes. 


1  TrTr'T;'i?r'''T'c  ov  a  TC''T'i7‘T>Rr4r4'v  cir  a  oir  r^ivT  a  irorirkA/NT  a  A/TTiT'c 

V*  AW  A  A  ^  A  A  4p~AA  a  a  <  Aw?  A  kwf  A  A  AW  'w' A  '«  A^.WA  VSW  A  A  ^ 


The  basic  results  of  this  project  are  contained  in  Figures  14  through  17.  These  Figures 
show  the  experimental  data  for  drag  coefficient,  slope  of  the  lift  and  moment  coefficient 
and  the  center-of-pressure  location  obteiined  from  the  range  tests.  The  Figures  also  show 
the  computed  results  from  the  design  codes  Euid  from  the  PNS  code. 

Figure  14  shows  the  zero-lift  drag  coefficient  for  all  eight  configurations  at  a  Mach 
number  of  4.0  and  a  Reynolds  number  based  on  model  cylinder  diameter  of  2.36x10'^.  The 
baseline  configuration,  (-1),  has  a  measured  drag  coefficient  of  0.56.  The  lowest  drag  is 
obtained  from  the  boattailed  configuration,  (-5).  This  is  not  surprising  since  the  boattail 
affects  the  base  pressure,  which  in  turn  determines  the  major  component  of  the  drag.  The 
largest  drag  is  produced,  as  expected,  by  the  biconic  shape  with  its  12  degree  flare  and 
consequent  high  wave  drag  component. 

The  prediction  techniques  determines  the  drag  coefficients  quite  well.  The  most  prob 
able  error  betw'een  the  NSWC  Aeroprediction  code  and  the  data  is  0.027  in  drag  coefficient. 
The  corresponding  value  for  the  PNS  results  is  0.029. 

Figure  15  shows  the  effect  of  the  different  configurations  on  the  pitching  moment 
about  the  center  of  gravity  of  the  projectiles.  The  center  of  gravity  varied  only  slightly 
between  configurations;  the  differences  in  moment  are  mainly,  but  not  entirely,  due  to 
aerodynamic  effects.  Obviously,  the  two  most  stable  configurations,  (-2)  and  (-4),  are  also 
the  configurations  with  the  highest  drag.  Configuration  (-6)  has  somewhat  higher  drag 
than  the  baseline,  but  it  has  much  improved  stability. 
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The  predictive  methods  genereilly  give  a  smaller  negative  pitching  momeiit.  The 
Second  Order  Shock  Expansion  method  consistently  gives  slightly  better  agreement  with 
the  experimental  data  than  the  PNS  computation. 


Figure  16  shows  the  lift  coefficient  slope  and  Figure  17  shows  the  center-of-pressure 
mee^sured  from  the  blunted  nose  of  the  projectiles.  These  two  figures  provide  an  alternative 
way  of  looking  at  the  pitching  moment  slope  results.  In  general,  the  Missile  DATCOM  code 
adequately  predicts  both  the  hft  slope  eind  the  center-of-pressure,  whereas  the  PNS  code 
consistently  underpredicts  lift  and  center-of-pressure.  It  is  possible  that  the  DATCOM 
results  are  fortuitous,  however,  because  the  DATCOM  code  predicts  the  forces  neglecting- 
viscous  effects.  The  smooth  body  viscous  effects  are  accounted  for  in  the  PNS  code  and 
the  afterbody  lift  is  significantly  less  than  predicted  by  the  design  code,  resulting  in  the 
more  forward  center-of-pressure.  However,  the  buttress  grooves  are  not  accounted  for  in 
the  computations.  One  possible  explanation  is  that  the  enlarged  boundary  layer  produced 
by  the  grooves  is  greatly  distorted  circumferentially  at  angle  of  attack  so  that  the  lee-side 
fljire  becomes  less  effective  relative  to  the  wind-side.  An  alternate  explanation  for  the 
lower  afterbody  lift  in  the  computations,  relative  to  the  experiment,  is  that  the  cavity  in 
the  rear  of  the  models  is  so  large  (over  4.0  calibers)  that  a  small  internal  pressure  difference 
between  the  top  and  bottom  surfaces  could  produce  a  considerable  lift  with  a  sufficiently 
large  moment  arm  to  account  for  the  observed  discrepancies.  CFD  studies  and  solid  base 
experiments  are  under  way  to  investigate  these  possibilities. 


Although  the  DATCOM  code  matches  the  data  fairly  accurately,  there  are  some 
discrepancies  between  predicted  and  observed  data  tieiids.  As  cjiu  be  seen  in  Figure  15, 
both  the  range  data  and  the  PNS  code  show  the  configuration  CS-V4-5  (boattail)  is  less 
stable  than  configuration  CS-V4-1  (baseline).  The  DATCOM  results  show  CS-V4-5  to 
be  more  stable.  This  trend  also  exists  for  the  lift  coefficient  slope  and  center-of-pressure 
location,  although  it  is  not  as  clearly  visible.  This  reversal  of  the  data  trends  can  lead  to 
serious  problems  when  choosing  a  projectile  geometry  for  a  particular  application,  which 
is  one  function  of  a  projectile  design  code. 


4.  PRESSURE  DISTRIBUTIONS 

One  of  the  primary  advantages  of  using  the  PNS  code  is  the  detail  which  can  be 
obtained  concerning  the  pressure  and  force  distribution  for  the  various  shapes,  from  which, 
the  overall  results  of  the  experiment  can  be  better  eissessed.  Figure  18  shov/s  the  wall 
static  pressure  for  the  CS-V2  configuration  at  three  Mach  numbers,  2.0,  3-0,  and  4.0.  The 
pressures  on  the  conical  nose  are  in  very  good  agreement  with  cone  flow  calculations.  The 
pressures  on  the  flares  are  also  approaching  cone  flow.  At  the  cylinder-flare  junction  the 
pressures  over-shoot  corresponding  to  local  two-dimensional  flow  conditions.  This  over¬ 
shoot  is  not  abrupt  because  of  the  fairing  used  to  prevent  oscillations  in  the  PNS  code. 
Figure  19  shows  the  pressure  distribution  on  the  leeward  side  of  the  first  six  configurations 
at  one  degree  angle  of  attack.  Among  the  first  five  configurations,  the  eftectiveness  of 
decreasing  the  slope  of  the  body  on  the  last  caliber  is  clearly  evident.  Configuration  (-6) 
produces  its  higher  stabilizing  lift  because  of  the  larger  flare  angle  (9.53°  compared  to  the 
6  degrees  of  the  others)  and  because  of  the  transition  to  a  square  cross  section  the  pressure 
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drops  to  nearly  the  same  level  as  the  cylindrical  skirt  body. 


5.  SKIN  FRICTION  DISTRIBUTION 

The  skin  friction  distribution  for  the  CS-V2  configuration  at  three  Mach  numbers  is 
shown  in  Figure  20.  At  a  Mach  number  of  2.0  there  is  a  strong  dip  at  the  cylinder-flare 
junction,  despite  the  0.25  caliber  fairing  at  that  point.  This  problem  almost  disappears  at 
the  higher  Mach  numbers.  The  flare  angle  of  four  degrees  is  so  slight  that  there  is  only  a 
small  increase  in  friction  on  the  aft  flare.  There  is  a  strong  Mach  number  effect  on  the  zero 
pressure  gradient  boundary  layer  on  the  cylinder  and  the  flaire  which  produces  a  decreasing 
skin  friction  drag  with  increasing  Mach  number. 

Figure  21  shows  the  effect  of  five  flare  configurations  on  the  skin  friction  distribution. 
This  Figure  is  for  Mach  number  of  4.0  and  zero  angle  of  attack.  The  PNS  code  predicts 
the  same  skin  friction  on  these  five  bodies  except  on  the  flare  extension.  The  6  degree  flare 
results  in  twice  the  increase  in  friction  as  that  obtained  on  the  4  degree  flare  of  the  CS-V2 
case.  The  12  degree  biconic  flare  has  the  greatest  friction  increase  whereas  the  friction  on 
the  boattail  extension  drops  below  the  cylinder  friction. 


6.  NORMAL  FORCE  DISTRIBUTION 

Figure  22  shows  the  distribution  of  the  normal  force  in  terms  of  the  slope  of  the  normal 
force  as  a  function  of  X/D  for  the  five  configurations.  The  nose  is  found  to  contribute  1.97 
to  the  slope  which  is  97  per  cent  of  the  slender  body  prediction.  The  cylinderical  mid  body 
contributes  only  about  half  of  that  of  the  nose  at  this  Mach  number  but  tends  to  provide 
an  increasing  contribution  with  Mach  number.  The  8.5  caliber  cylinder  is  sufficiently  long 
so  that  increeising  its  length  would  not  provide  significant  additional  lift.  The  6  degree,  3.5 
caliber  flare  contributes  more  lift  than  the  entire  forebody.  The  12°  biconic  substantially 
increases  the  normal  force  and  stabilizing  moment  but  at  the  cost  of  an  increase  in  drag. 
The  boattailed  extension  (-5)  produces  only  a  little  less  added  lift  than  the  cylindrical 
extension  (-3)  but  with  an  important  decrease  in  base  drag. 

Figure  23  shows  a  comparison  between  the  normal  force  predicted  by  the  Shock  Ex¬ 
pansion  method,  used  in  the  DATCOM  and  NSWC  codes,  and  the  PNS  code.  It  shows 
that  the  consistently  higher  prediction  by  the  SOSE  method  is  because  it  predicts  more 
lift  on  the  flare.  That  method  is  an  inviscid  one  whereas  the  PNS  result  takes  the  viscous 
flow  into  account  with  a  reduced  efficiency  of  the  flare  (85  per  cent  less  normal  force  for 
the  4.5  caliber,  6  degree  flare).  This  is  further  confirmed  in  comparing  the  circumferential 
pressure  distribution  as  shown  in  Figure  24.  The  pressure  distributions  are  for  a  point 
half  the  flare  length  and  the  predicted  pressures  levels  are  significantly  different.  sub¬ 
tracting  the  pressure  coefficients  at  90°,  the  difference  in  the  windward  and  leeward  sides 
which  determines  the  normal  force,  is  displayed.  Although  the  SOSE  method  gives  a  larger 
difference,  the  distributions  are  similar. 
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7.  EVALUATION  OI  THE  MOST  EFFECTIVE  CONFIGURATION 


Figure  25  presents  the  meEisured  data  as  a  correlation  of  drag  coefficient  with  center 
of  pressure  location.  This  format  is  chosen  to  emphasize  the  trade  offs  between  stability 
and  drag.  The  ’’best”  configurations  have  the  most  rearward  center  of  pressure  and  the 
least  drag,  and  are  thus  located  toward  the  bottom  right  corner  of  the  figure.  A  single 
continuous  curve  can  be  drawn  through  the  data  of  configurations  (-1)  and  (-2)  showing 
the  effect  of  flare  length  for  a  single  flare  angle,  i.e.,  increasing  the  flare  length  from  3.5 
calibers  (configuration  (-1))  to  4.5  calibers  (configuration  (-2))  increases  the  drag  coefficient 
by  0.12  and  shifts  the  center  of  pressure  back  by  1.1  erdibers.  Configurations  (-5),  (-3), 
(-2)  and  (-4)  also  form  a  related  series  of  afterbody  shapes,  starting  with  the  boattail  and 
ending  with  the  12  degree  biconic  body.  A  curve  drawn  through  these  four  points  shows 
the  effect  of  the  angle  of  a  one  caliber  extension  to  the  baseline  flare.  The  figure  shows 
that  flare  extensions  with  larger  angles  than  the  baseline  flare  are  less  efficient  than  simply 
extending  the  flare.  Conversely,  extensions  with  a  shallower  flaire  angle  axe  desirable.  The 
fin  configuration  (-7),  is  judged  the  most  effective  and  (-6)  and  (-8)  are  the  next  most 
efficient  shapes  for  improving  stability  at  a  minimum  penalty  in  drag. 


V.  CONCLUSIONS 

Ten  flare-stabilized  projectile  configurations  have  been  tested  in  the  BRL  Aerody¬ 
namics  Ballistic  Remge  and  the  results  compared  to  numerical  computations  based  on  two 
design  codes  and  a  Parabolized  Navier-Stokes  code.  The*  evaluation  of  the  results  show: 

The  configuration  with  a  combination  of  6  degree  flare-boattail  with  small  fins  pro¬ 
duced  the  most  stable  shape  with  the  least  drag  penalty  as  compared  to  conventional 
flares. 

The  prediction  of  the  aerodynamic  characteristics  using  the  NSWC  Aeroprediction 
code  gave,  in  general,  very  good  results,  despite  the  fact  that  effective  bodies  had  to  be 
devised  to  represent  the  actual  body  shapes,  and  that  corrections  had  to  be  introduced  to 
account  for  the  buttress  thread  roughness  effects.  However,  some  of  the  data  trends  for 
biconic  configurations  were  not  predicted  correctly. 

The  Parabolized  Navier-Stokes  results  complemented  the  sem.i-empirical  computations 
by  providing  more  detail  on  the  development  of  the  aerodynamic  forces.  The  PNS  pitching 
moment  prediction  is  consistently  less  than  the  measurements,  possibly  because  of  the  effect 
of  the  open  base  cavities  or  of  the  buttress  groove  boundary  layer  on  the  projectiles  tested. 


VI.  RECOMMENDATIONS 


1.  It  is  recommended  not  to  use  the  Missile  DATCOM  code  for  the  evaluation  of  the 
axial-force  coefficient  of  bodies  with  flares,  since  it  over-predicts  the  contnbution  of  the 
base. 
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2.  It  is  recommended  to  use  the  Missile  DATCuM  code  for  predicting  the  normal- 
force  and  the  pitching-moment  coetHcients  of  bodies  in  the  transonic  regime.  The  reason 
being  the  accuracy  of  its  aerodynamic  characteristics  data  base  in  this  Mach  number  range. 

3.  It  is  recommended  to  use  the  option  of  truncated  nose,  in  the  DATCOM  and 
NSWC  codes,  only  for  configurations  having  non-small  ratio  of  truncation  diameter  to 
maximum  diameter.  Cases  having  a  slightly  truncated  nose  should  be  emalyzed  as  if  the 
nose  is  extended  to  a  point. 

4.  It  is  not  recommended  to  use  the  NSWC  or  DATCOM  codes  for  design  opti¬ 
mization  of  biconic  flared  projectiles.  These  codes  do  not  correctly  predict  the  stabiliz- 
ing/destabiiizing  effect  of  the  second  flare  section. 
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TABLE  2:  TYPICAL  TABULATED  XM910  RESULTS 


MACH  NUMBER=2.0 

■Si 

XcpID 

CS-V2 

CS-V4-1 

CS-V4-2 

CS-V4-3 

CS-V4-4 

CS-V4-5 

0.150 

0.227 

0.270 

0.224 

0.528 

0.303 

0.191 

0.180 

0.201 

0.200 

0.202 

0.197 

0.414 

0.395 

0.496 

0.336 

0.698 

0.211 

0.090 

0.090 

0.090 

0.090 

0.090 

0.090 

0.835 

0.892 

1.057 

0.850 

1.518 

0.801 

5.58 

5.95 

7.52 

6.70 

8.56 

6.11 

9.26 

9.29 

10.43 

10.07 

11.45 

9.47 

MACH  NUMBER=3.0 

CS-V2 

CS-V4-1 

CS-V4-2 

CS-V4-3 

CS-V4-4 

CS-V4-5 

CS-V4-6 

0.122 

0.180 

0.214 

0.17S 

0.408 

0.217 

0.273 

0.323 

0.406 

0.297 

0.539 

0.208 

0.080 

0.080 

0.080 

0.080 

0.080 

0.080 

0.080 

0.646 

0.746 

0.882 

0.733 

1.213 

0.680 

5.75 

5.76 
7.15 
6.22 
8.49 
5.70 
7.98 

7.76 

8.64 

10.12 

9.20 

11.12 

8.52 

10.39 

MACH  NUMBER=4 

CS-V2 

CS-V4-1 

CS-V4-2 

CS-V4-3 

CS-V4-4 

CS-V4-5 

CS-V4-6 

0.105 

0.155 

O.ISO 

0.154 

0.348 

0.177 

0.155 

0.148 

0.167 

0.162 

0.174 

0.158 

0.187 

0.217 

0.273 

0.202 

0.358 

0.144 

J 

0.060 

0.060 

0.060 

0.060 

0.060 

0.060 

0.060 

0.507 

0.580 

0.686 

0.578 

0.940 

0.539 

5.47 

5.79 

7.03 

6.08 

8.68 

5.63 

7.43 

8.13 

8.33 

9.73 

8.70 
11.01 

8.17 

9.71 

MACH  NUMBER=4.6 

CSV2 

0.100 

0.110 

0.150 

0.050 

0.410 

5.94 

7.94 

14 
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7.8r 
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3.27 


1.00 

-I,. 


8.00* 


MI 


T 


M 


0.45 


9.98 


11.80 


T 


15.36 


ALL  DIMENSIONS  IN  CALIBFiS  (1  CALIBER  =  8. 28  mra ) 


Figure  1.  Sketch  of  baseline  projectile  with  CS-V4-1  flare  afterbody. 
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Figure  2.  Afterbody  and  extension  configuration. 
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Figure  4.  Spark  shadowgraph  of  CS-V4-1  baseline  projectile  in  flight  (high 
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Figure  5.  CoiiTpaxison  between  predicted  characteristics  of  a  cone-cylinder  bo<i} 
and  Barth’s  data  for  the  slope  of  normal-force  curve. 
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Figure  6.  Comparison  between  predicted  characteristics  of  a  cone-cylindei  body 
and  Barth’s  data  for  the  center-of-pressure  location. 
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Figure  8.  Comparison  between  computed  forebody  charactoristics  and  eiu]niica.l  data 
for  the  location  of  center-of-pressure. 
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Figure  9.  Geometry  of  equivalent  flare  stabilizers. 
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Figure  10.  Axial-forre  coeffiriont  a.s  a  fimction  of  Macli  miiriber 
f  -  the  CS-Vl  configuration. 


Figure  11. 

Axial-force  coefficient  as  a  fiinction  of  Mach  number 
with  component  brakedown  lor  the  CS-V2  confiffuration. 
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Figure  13.  Comparison  between  analysis  and  ranse  results  for  the  location  of  the 
center-of-pressure  as  a  function  of  Mach  number. 
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Figure  14.  Zeio  lift  drag  coefficient  for  all  confip;urations  at  Moo~4.0. 
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Figure  15.  Slope  of  the  pitching;  moment  coefficient  about  the  center  of  )j;ravity 
at  Moo=4.Q. 
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Figure  17.  Center  of  pressure  location  from  the  projectile  nuse  for  all  configurations 
at  Mo»=4.Q. 
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Figure  21.  Skin  friction  coefficient  distribution  for  five  configurations  at  M. 


Figure  22.  Normal-force  distribution  for  six  confiRurations  at  A4o=4.0  g  =  0^. 
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Figure  23.  Comparison  between  Second  Order  Shock  Expansion  and  PM.q 
■force  distgihiit.inp,,  - - - 
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Figure  25.  Correlation  of  drag  coefficient  with  center  of  i:)ressurR  location  at  Mqq -4.0. 
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LIST  OF  SYMBOLS 


A 

lateral  surface  area 

b 

span  of  the  fin  or  straJce 

Cd 

drag  coefficient 

Cl 

skin  friction  coefficient 

C, 

pressTore  coefficient 

Cl 

lift  coefficient 

Cn 

normal-force  coefficient 

A  Cn 

contribution  of  flaxe  to  normal-force  coefficient 

Cm 

pitching  moment  coefficient 

D 

projectile  diameter 

Dk 

base  diameter 

d 

width  of  a  side  of  a  square  body 

M 

Mach  number 

S 

body  local  cross  sectional  area 

Sbt 

surface  area  covered  with  buttress  threads /  grooves 

Sr 

reference  area, 

X 

distemce  from  projectile  nose  along  the  axis  of  symmetry 

Xcp 

location  of  the  center  of  pressure 

Greek  Symbols 

a 

Angle  of  attack 

rj 

flare  efficiency 

P 

density 

4> 

circumferential  angle 

<t>j 

flare  angle 

Subscripts 

b 

base 

bg 

buttress  groove 

e 

equivalent  slender  body 

f 

viscous  friction 

0 

zero  lift  conditions 

oo 

conditions  at  infinity 
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